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A rigorous derivation of quantum Langevin equation from microscopic dynamics in the 
low density limit is given. We consider a quantum model of a microscopic system (test 
particle) coupled with a reservoir (gas of light Bose particles) via interaction of scattering 
type. We formulate a mathematical procedure (the so-called stochastic golden rule) 
which allows us to determine the quantum Langevin equation in the limit of large time 
and small density of particles of the reservoir. The quantum Langevin equation describes 
not only dynamics of the system but also the reservoir. We show that the generator of 
the corresponding master equation has the Lindblad form of most general generators of 
completely positive semigroups. 
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1. Introduction 

In the past years many important physical models have been investigated by using 
the stochastic limit method (see [1, 2] for a survey). In particular Fermi golden 
rule has been generalized to a stochastic golden rule which allows to solve the fol- 
lowing problem: give a quantum Hamiltonian system, determine the associated 
Langevin equation (including the structure of the so-called master fields or quan- 
tum noises) and obtain as a corollary the master equation associated to this sys- 
tem. 

The Langevin equation is an asymptotic equation, valid in a regime of large 
times and small parameters (e.f. weak coupling) and its interest for the applications 
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comes from the fact that the relevant physical phenomena in this regime are much 
more transparent than in the original Hamiltonian equation and in this sense an 
approximate equation gives a better description of the physical phenomenon than 
the original exact equation. 

The stochastic golden rule plays a fundamental role in the stochastic limit be- 
cause it allows to separate the problem of finding the correct solution from the 
problem of giving a rigorous proof of the convergence of the method. 

Although well developed in the weak coupling case, in the low density case, the 
stochastic golden rule was formulated in [1] only for Fock fields and this limitation 
exclude of the most interesting cases for the applications. 

The development of a stochastic golden rule for the low density limit required a 
totally different approach from that developed for the weak coupling case and the 
first step in this direction was done in the paper [3] where the problem was solved 
under the (very strong) assumption of rotating wave approximation (RWA). In the 
present paper we bring to a completion the programme begun in [3] and we establish 
the stochastic golden rule for general discrete spectrum system interacting with a 
Boson reservoir without RWA. This also gives, when combined with the analytical 
estimates of part III of [1], a simplified proof of the deduction of the stochastic 
Schrodinger equation for the low density limit. The new proof is much simpler and 
allows a deeper understanding of the emergence of the scattering operator, more 
precisely the T-operator, in this equation. 

To describe a quantum physical model to which we apply the low density limit let 
us consider an N-level atom immersed in a free gas whose molecules can collide with 
the atom; the gas is supposed to be very dilute. Then the reduced time evolution for 
the atom will be Markovian, since the characteristic time ts for appreciable action 
of the surroundings on the atom (time between collisions) is much larger than the 
characteristic time Ir for relaxation of correlations in the surroundings. Rigorous 
results substantiating this idea have been obtained in [4] . 

The dynamics of the N-level atom interacting with the free gas converges, in 
the low density limit, to the solution of a quantum stochastic differential equa- 
tion driven by quantum Poisson noise [5]. Indeed, from a semiclassical point of 
view, collision times, being times of occurrence of rare events, will tend to become 
Poisson distributed, whereas the effect of each collision will be described by the 
(quantum-mechanical) scattering operator of the atom with one gas particle (see 
the description of the quantum Poisson process in [6]). 

We consider a microscopic system (test particle) interacting with a gas of Bose 
particles (reservoir) via interaction of scattering type so that the particles of the 
reservoir are only scattered and not created or destroyed. The reservoir is supposed 
to be a very dilute gas at equilibrium state. We obtain starting from microphysical 
dynamics, in the limit of long time and small density of Bose particles, the quantum 
Langevin equation which describes the evolution of any system observable. The 
Langevin equation includes not only the system but also the reservoir dynamics. 
After averaging over the reservoir equilibrium state we obtain a master equation 
which describes the evolution of any system observable, i.e. the reduced evolution of 
the test particle. The generator of this master equation has the standard Lindblad 
form [7]. 
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Starting from the quantum stochastic equation (7.29) for the cvohition operator 
one can easy deduce the corresponding quantum hnear Boltzmann equation for the 
density matrix of the system. This is a Hnear dissipative equation describing an 
irreversible evolution of the test particle. A general structure of a master equa- 
tion driving a completely positive time evolution for a test particle in a quantum 
(Rayleigh) gas has been investigated in [8]. It turns out that general requirements 
of the underlying symmetry arc not sufficient to fix the form of the master equation 
and one needs a microscopic derivation of the master equation [9]. A linear dissi- 
pative equation describing the long time dynamics of a particle interacting with a 
reservoir (a quantum Fokkcr-Planck equation) has been studied in [10] where the 
problem of rigorous derivation of such equations was emphasized. 

Using the white noise approach developed in [1] and the energy representation [3] 
we derive a new algebra for the quadratic master fields in the low density limit. An 
advantage of this method is the simplicity of derivation of the white noise equations 
and in the computation of correlation functions. We show that the drift term in the 
QSDE is simply related with T-operator describing the scattering of the system on 
the one particle sector of the reservoir. 

The main results of the paper are: 

(i) the explicit determination of the algebra of the master field in energy repre- 
sentation (section 3) 

(ii) the determination of the white noise Hamiltonian equation (6.20) for the low 

density limit 

(iii) the determination of the (causally) ordered form (7.29) of this equation 
which, as known from the general theory developed in [1], is equivalent to a quantum 
stochastic differential equation in the sense of [11] 

(iv) formula (8.37), which gives an explicit expression for the drift in terms of 
the 1-particle T-operator 

(v) the determination of the quantum Langcvin equation (9.39) and the quantum 
Markovian generator (9.40) of the corresponding master equation. 

Our strategy is the following. In section 2 we describe the model which we 
consider; and for this model, we introduce the Fock-anti-Fock representation. In 
section 3 we introduce the rescalled fields whose stochastic limits are the master 
fields in the LDL. In section 4 we derive the commutation relations of the master 
fields (Theorem 1). Then in section 5 we construct a concrete representation of this 
algebra. Using this representation we write in section 6 the stochastic Schrodinger 
equation (6.20) for the limiting evolution operator. In section 7 we rewrite the 
stochastic Schrodinger equation in normally ordered form, i.e. in the form of a 
quantum stochastic differential equation. In section 8 we prove formula (8.37) which 
expresses the drift, in the equation for the evolution operator in the LDL in terms 
of: the 1-particle T-operator, the 1-particlc gas Hamiltonian and the spectral pro- 
jections of the system Hamiltonian. In section 9 we obtain the quantum Langevin 
equation and the quantum Markovian generator of the corresponding master equa- 
tion. 

2. The Model 

Now let us explain our notations. Let Hs be a Hilbert space of the system 
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with the Hamiltonian H^. We will suppose that the system Hamiltonian has a 
discrete spectrum 

k 

where ek are the eigenvalues and Pk the spectral projections. 

We will consider the case of a Boson reservoir in this paper. Therefore the 
reservoir is described by the Boson Fock space r(7ii) over the one particle Hilbert 
space Til = L^(R'^) (the scalar product in Tii we denote as < •, • >), where d = 3 
in physical case. Moreover, the free Hamiltonian of the reservoir is given by := 
r(iJi) (the second quantization of the one particle Hamiltonian Hi) and the total 
Hamiltonian of the compound system is given by a self-adjoint operator on the total 
Hilbert space Hs T{Hi), which has the form 

Htot := Hs®l + l® Hb + V =: Ho + V 

Here V is the interaction Hamiltonian between the system and the reservoir. The 
evolution operator at time t is given by: 

U{t) := ej*"° ■ e-'*"'°' (2.1) 

and it satisfies the differential equation 

dtU{t) = -iv{t)u{t) 

where the quantity the evolved interaction V{t) is defined by 

V{t) = e'*"°Ve-'*"° . 
The interaction Hamiltonian will be assumed to have the following form: 
y = £) A+{go)A{gi) + D+ (g) A+{gi)A{go) 

where is a bounded operator in Hs, D G B(Hs), A and are annihilation and 

creation operators and go,9i G Hi are form- factors describing the interaction of the 
system with the reservoir. Therefore, with the notations 

St := e''"' ; D{t) := e''"' De''*"' 

the evolved interaction can be written in the form 

V{t) = D{t) ® A+{Stgo)A{Stg{) + D+{t) A+{Stgi)A{Stgo) (2.2) 

The initial state of the compound system is supposed to be of the form 

P = PS® <fi^^^- 

Here ps is arbitrary density matrix of the system and the initial state of the reservoir 
</j^^^ is the Gibbs state, at inverse temperature /?, of the free evolution, i.e. the gauge 
invariant quasi-free state characterized by 



^(«)(W(/)) = exp(-i </,(! + ee-^«^)(l - ee-^^^)-V >) (2-3) 
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for each / G Jii. Here W{f) is the Weyl operator, (3 the inverse temperature of 
the reservoir, ^ = e^^^ the fugacity, /U the chemical potential. We suppose that 
the temperature (3~^ > 0. Therefore for sufficiently low density one is above the 
transition temperature, and no condensate is present. 

We will study the dynamics, generated by the Hamiltonian (2.2) and the initial 
state of the reservoir (2.3) in the low density regime: n — > 0, t ~ 1/n (n is the 
density of particles of the reservoir) . We do not fix the initial state of the system 
so our results can be applied to an arbitrary initial density matrix of the system. 
In the low density limit the fugacity ^ and the density of particles of the reservoir 
n have the same asymptotic, i.e. 

limiM = i 

n^O n 

Therefore the limit n ^ is equivalent to the limit ^ ^ 0. 

Throughout this paper, for simplicity, the following technical condition is as- 
sumed: the two test functions in the interaction Hamiltonian have disjoint supports 
in the energy representation. This condition is invariant under the action of any 
function of H\ and implies that the two test function go, gi in the interaction Hamil- 
tonian satisfy: 

< go:Ste-^^"^gi>=0 € R. 

In the paper [3] we obtained for the model described above a quantum stochastic 
differential equation under additional rotating wave approximation condition: 

D{t) = e**^^L'e-'*'^^ = e^'^'^^D (2.4) 

for some real number luq. 

In the present article we will derive the white noise Schrodinger equation without 
assuming any relation between D and Hg. 

Let us rewrite the free evolution of in a form which is convenient for derivation 
of the white noise equation. For this we introduce the set of all Bohr frequencies 

B(Hs) =: B := {w I 3 £fe, e„ e spec Hs s.t. u} = Sk - £m} 

Using this notion and the properties of the spectral projections one can rewrite the 
free evolution of an arbitrary system operator D in the form 

D{t) =J2^''^"''~""'^PkDPm = J2 e"""^a, 

where we denote 

:= Yl PkDPm 

k.7n:€m—ek—^ 

We realize the representation space as the tensor product of a Fock and anti- 
Fock representations. Then the expectation values with respect to the state tp^^' for 
the model with the interaction Hamiltonian (2.2) can be conveniently represented as 
the vacuum expectation values in the Fock-anti-Fock representation for the modified 
Hamiltonian. 
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Denote by T-fi the conjugate of Hi, i.e. is identified to H\ as a set and the 
identity operator l : Hi — Hi is antihnear 

t(A/) := A.(/) 

<L{f),L{g) >,~<g,f> 

Then, Wi is a Hilbcrt space and, if the vectors of Hi are thought as ket-vectors 
1^ >, then the vectors oiH'i can be thought as bra-vectors < ^|. The corresponding 
Fock space T{H[) is called the anti-Fock space. 

It was shown in [5] that, with the notations Do^.o = D,^, D^^i = Z)+ and A = 
the part of the modified Hamiltonian which gives a nontrivial contribution in 
the LDL, acts in Hs O r(Wi) T{H'i) as 

^A(i) = E E ^-,.e-(-i)^'*- ® (A+{Stg,)A{St9i-e) ® 1 

e=0,l weS 

+A(A(5t5i-.) ® A{StLg,) + A+{Stge) ® A+ {StLgi_,))) . 

Here A and are Bose annihilation and creation operators acting in the Fock 
spaces T{Hi) and r(Wi) and L := e''^"^/'^. 

The interaction Hamiltonian Hx (t) determines the evolution operator t//'^^ which 
is the solution of the Schrodingcr equation in interaction representation: 

with initial condition 

= 1. 

This is equivalent to the following integral equation for the evolution operator 

t 



The convenience of the Fock-anti-Fock representation, as it was mentioned above, 
is in the fact that the expectation value, for example, of any time evolved system 
observable Xt = U*{t){X (g) l)U{t) with respect to the state (2.3) is equal to the 
vacuum expectation value in the Fock-anti-Fock representation, i.e. 

ip^^\U*{t){X ^ l)U{t)) =< U;^^\X ® 1 ^ l)C/f ^ >vac 

3. Energy Representation 

We will investigate the limit of the evolution operator uj:^^ when ^ +0 after the 
time rescaling t t/^, where ^ = A^. After this time rescaling the equation for the 
evolution operator becomes 

9tu^/l = E E ® [Ne,i-eAi-'^y^,t) 

weBe=0,l 
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+Bl_,,e,A((-l)'-^a;,i) + Bj;i-e,A((-ir'^.*)]f^t^/i 

where we introduced for each ei, £2 = 0, 1 and w G B the rescaled fields: 

N,„,,,^{uj,t) := ^e-''"/^'A+{St/;,2g,,)A{St/x-9e,) ® 1 (3.5) 



and -B^ xi^^t) is the adjoint of Bs^^s2,\{'^tt)- 

The energy representation for the creation and annihilation operators is defined 
by the formulae 

A+{g) := A+{PEg), Ae{9) := A{PEg) (3.7) 

with 

Pe:=^ J dtSte-'*'' = S{Hi - E) (3.7a) 
—00 

It has the properties 

dEPse'*^ (3.76) 

PePe' = S{E - E')Pe, P*e = Pe 

The meaning of the (5-function in (3.7a) is explained in [1], section (1.2). In our 
case for Tii = L'^{W^) the one-particle Hamiltonian is the multiplication operator 
by the function uj{k) and acts on an element / e Li^iW^) as {Hif){k) = u){k)f{k) 
so that Pe = 6{uj{k) - E). 

It is easy to check that 

[AeU). A%,{g)] = 5{E - E') < f, Psg > . 
Using the energy representation (3.7b), (3.5) and (3.6) become respectively 

Ar,,,,„A(a;,t) = j J dEidE2N,,^,,,xiEi,E2,uj,t) 

B,„,„A(w,t)= j j dEidE2B,,,,,^x{Ei,E2,oj,t) 
where the rescaled fields in energy representation are given by 

git(Ei-£;2-u;)/A2 

N,,,,,^x{El,E2,iJ,t) := A+^{9er)AEA9e,) 1 (3.8) 

^it(E2-Ei+uj)/X^ 

B,,^,,^x{Ei,E2,uj,t) := AE,{ger)<^AE,{Lg,,). (3.9) 
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and -B^ ^{Ei, E2,u),t) is the adjoint of B^^^s:2,\{E\,E2-,uJ,t). Let us also denote 



^,{E) := j <g,,Stg,>e-''^dt. 

— oo 

4. The Algebra of the Master Fields 

In this section we derive the algebra of master fields which is the limit as A —> of 

rescalcd fields (3.8), (3.9) in the sense of convergence of correlators. 

Denote Vt the lattice of all linear combinations of the Bohr frequencies with 
integer coefficients: 

N 

O = {w I w = rifeW/c with A/' e N, rife e Z, Wfe e B} 

k=l 

and extend the definitions (3.8) and (3.9), of the rescaled fields, by allowing in them 
an arbitrary a; e O. The following theorem describes the algebra of commutation 
relations for the master field in the LDL. 

Theorem 4.1 The limits of the rescaled fields 

Xe, ,e, {Ei,E2,u;,t):= \im X,,,,,,x{Eu E2,u;,t), X = B,B+,N 

exist in the sense of convergence of correlators and satisfy the commutation relations 
[B,,,,,{EuE2,iJ,t),B+^^^{E3,Ei,uj',t')] = 2Tr6^,^,S,,,,J,,,,J{t' - t) 

X 5{Ei - Es)5{E2 - Ei)5{Er -E2-uj)< g,„PE,ge, >< ge^Ps^L'^g,^ > (4.10) 

[Be„e,{EuE2,lV,t),Ne„eAE3,E4,lv',t')]= 2TTde, ,eAt' " t) 

X 5{Ei - E3)5{Ei -E2-w)< g,„PE,9e, > B,,^,,{Ei,E2,L0 - w' ,t') (4.11) 

[Af,„e,(Ei, E2, w, t), N,,,,,{E3, £4, w', t')] = 2nS{t' - t) 
X i^5e^^eME2 - E2.)5{E:i - El + uj) < g^^,PE^ge., > N,,^,^{Ei,E4,u+u',t') 

-5e„eAEi-Ei)5{Ei-Ei-J)<g^„PE,ge, >iVs3,e2(i^3,i^2,w + a;',t)} (4.12) 

The causal commutation relations of the master field are obtained replacing in 
(4-10), (4-11), (4-12) the factor 5{t' — t) by S+{t' — t) where the causal 5-function 
6+{t' —t) is defined in section (8.4) of [1], 2wd{Ei — E2—oj) by {i{Ei — E2—oj — iO))~^ 
and 2-k5{E2. - iJi ± w) by {i{E3 -Ei±ij- iO))-^ 
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Proof. Introduce the operators: 

git(_E2-Bi-w)/A^ 

Nl^^^^^^{Ei,E2,io,t) := l®A%^{Lg,,)AEALge^) 

with AsiLgs) and A'^{Lgs) defined by (3.7), the commutators of the rescalled fields 
are: 

„i{t'-t){Ei-E2-iv)/\^ 

X (d,,,,J,,^,,e'*'^^-'-'^/^"S{E^ - E3)6{E2 - £4) < 9.,,PE,9e. >< 9e.,PE,L^9e. > 
+ X^6,,,,J{Ei - E3) < g,„PE,ge, > Nl^^^^^^{E^, E2,io' - Lo,t') 
+ \^5,,,,J{E2 - Ei) < g,„PE,L''g,, > N,,,,,^x{E3,Ei,uj' - w.t')) (4.13) 

^i{t'-t)(Ei-E2-uj)/X'^ 

[B£i,e2,\{El,E2,(J,t),Nss,e4,\{E3,E4,Ul',t')] = 

x5,,^,J{E,-E3) <g,„PE,ge,> B,,^,,^x{E4,E2,L0-u;',t') (4.14) 

[N,,,e2AEl,E2, W, t), JV,3,e4,A(^3, E4, J, t')] 
gi(t'-t)(B3-Si+u;)/A^ 



A2 

gi(('-/)li:;:i-i'l-uj' )/A- 



5e2,eME2 - E3) < 9e2,PE29e2 > N,,,e^,x{El, Ei,UJ + J ,t') 

^2 6e^,eJ{El - Ei) <ge^,PE^gei > Ne^^e^^Ez, E2,U1 + Uj\t) . 

(4.15) 

Notice that in the sense of distributions one has the limit 

hm '"^'o^' ' e^*'("-'^')/^' = 2T,5^,^,5{t' - t)5{E^ - E2 - uj) (4.16) 

and, in the sense of distributions over the standard simplex (see [1] for details) one 

has the Umit 

i{t'-t){E^-E2-uj)/\' 1 

lim e'*(--)/A ^5 s ^i>_t^ (4.17) 

A^o A^ i[Ei — E2 — oj — ivi) 

The proof of the theorem follows by induction from the commutation rela- 
tions (4.13)-(4.15) using the limits (4.16) and (4.17) and standard methods of the 
stochastic limit. □ 

5. The Master Space and the Associated White Noise 

In this section we construct a representation of the limiting algebra (4.10)-(4.12). 



A Stochastic Golden Rule for the LDL 



Let Khc a, vector space of finite rank operators acting on the one particle Hilbert 
space Til introduced in section 2, with the property that, for any w G f2 and for 
any X, Y & K 



< X,Y>^:= j dtTr {e'l^"^ X* StY S;) 
= 27ry" dETr{e-f^"^X*PEYPE-S) < 



Because of our assumptions on Hi, the space K is non empty and < ■, ■ >u, defines 
a prescalar product on K. We denote {K, < ■, ■ >„} or simply K^) the Hilbert space 
with inner product < •, • >uj obtained as completion of the quotient of K by the 
zero < •,• >a;-norm elements. Denoting Yt = StYSf the free evolution of Y, one 
can rewrite the inner product as 



<X,Y >^= j dtTr {e-f^"'X*Yt) 



With these notations the representation space of the algebra (4.10)-(4.12) is the 
Fock space 

T{L\R+) ® if^) = r((g) L\R+,K^)) = (g) T{L\R+,K^)) 

where the infinite tensor product is referred to the vacmim vectors. With these 
notations the master fields are realized as a family of white noise operators &t,aj(') 

which act on r(L^(IR+) (g) K^^) and satisfy the commutation relations 

wen 

[btAX),bl^,{Y)] = d{t' - t)6^,^, <X,Y>^. 

Moreover each white noise operator 6t,a,(-) acts as an usual annihilation operator in 
r(L^(R_i_) ig) K^) and as identity operator in other subspaces. 

We will construct a representation of the algebra (4.10)-(4.12) in the Fock space 
defined above by the identification of the operators B^^^AEi,E2,L0,t) with the 
white noise operators 

B^^,AEuE2,u:,t) = bt,u,{\PE,ge, >< PE^Qe^) 

The number operators will then have the form 

N,^,,^{Ei,E2,L0,t) = 

£e{o,i}wien 
with 

l^e{E) :- ^ 



<g„PEL^ge>' 

One easily checks that these operators satisfy the algebra (4.10)-(4.12). 
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Let us introduce, for simplicity of calculations, the operators 

B,,^,,{E,uj,t) := j dE'B,,,,,{E',E,uj,t) 

N,,,,,{u,t) := j dEidE2N,,,e,{EuE2,uj,t) = j dEN,,,e,{E,io,t) 

with 

,OJi - U},t) 

e6{0,l}aJien 

The operators Bsi^e^{E,uj,i), Bf^^^^{E,ui,t) and N^-^^s^i^Tt) satisfy the (causal) 
commutation relations: 

{B,,,,,{E,u,,t),B+^^^{E',u',t')] 



= 5+{t' - t)5e,,eA.,eM^'5{E - E')-f,, {E + u)p^-^{E) (5.18) 

[Bei,e.(S,a;,t),7Ve3,e,(a;',i')] = 5+{t'-t)5e,,e,leAE+'^)B,,,eAE,uj-J,t'). (5.19) 

In these notations the limiting white noise Hamiltonian acts on r(L^(IR.+) (g) 

K^) as 

wen 

H{t)=Y^D^(^\^No^^{w,t)+ I dE\Bi^o{E,-w,t)+B+^{E,iv,t)\^ 
+ Y,Dt^(^[N^fl{^^,t) + j dE[Bo,^{E, -a;,i) + B+o(^,a;,i)]}. 



6. White Noise Stochastic Schrodinger Equation 

The results of the preceding section allow us to write the white noise Schrodinger 
equation for the evolution operator in the stochastic limit 

dtUt = -iH{t)Ut = -iJ2{Du.^[^o,i{^^t)+ I dE(^Bifi{E, -w,t)+B^^i{E,oj,t)^ 



ujeB 



+Dt. 



,o{iv,t) + J dE(^Bo,i{E,-u,t) + B+o{E,iv,t}'^'\}Ut (6.20) 

Following the general theory of white noise equations, in order to give a precise 
meaning to this equation we will rewrite it in the causally normally ordered form in 
which all annihilators are on the right hand side of the evolution operator and all 
creators on the left hand side. After this procedure we obtain a quantum stochastic 
differential equation (QSDE) in the sense of Hudson and Parthasarathy [11]. 
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For each oj G r2\B define = 0. Then define for any w, w' e the operators 
TO^,(iJ) andr^_,,(i?) by 

Tt,AE) :=^ME + t^ho{E + uj")D+,_^D^„_^, 

uj" 

We will also denote for each e = 0,1 

The following lemma plays an important role in the derivation of the normally 
ordered white noise equation. 

Lemma 6.1 In the above notations one has 

^(1 + 7o)^,^' {E)Bo,o{E, to', t)Ut = UtBo,o{E, uj,t)-iJ2 D^-^no{E + u) 

X f/tBi,o(S, uj',t)-J2 D^,+^D+ix^\E)^o{E + uj)^i {E - J)Ut (6.21) 

Uj' 

5^(1 + To)^,^' {E)Bo,i {E, u', t)Ut = UtBo,i{E, - i ^ D^_^,-io{E + uj) 

X UtB^,i{E,uj' ,t) - iD^iji-\E)-^o{E + uj)Ut (6.22) 
^ (1 + Ti )„,„. (i;)Bi,i {E, J, t)Ut = UtBi,,{E, u>,t)-iJ2 {E + oj) 

xUtBo,i{E,uj',t) -J2Di,_^D^'Hi\E)^i{E + ij)^o{E + uj')Ut (6.23) 

uj' 

^(1 + Ti)^,u,' {E)B^,o{E, uj', t)Ut = UtB,,o{E, cj,t)-iY, D^'-u^ME + oj) 

X UtBo^oiE,tu',t) -iD+^Ho\E)jiiE + Lu)Ut (6.24) 

Proof. It is clear that 

B,, ^,,{E,u,t)Ut = {B,, {E,uj,t),Ut] + UiB,,,,,{E,ij,t) 

Using the integral form of equation (6.20) for the evolution operator one gets 

t 

[B,,,,,{E,Lo,t),Ut]=-i J dti([B,,^,,iE,u;,t),H{ti)]Ut, 
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+H{ti)[B,,^,,{E,LJ,t),Ut,]) (6.25) 

Notice that due to the time consecutive principle (see [1] for details) one has for 
t > ti 

Using the causal commutation relations (5.18), (5.19) one can compute the 
causal commutator [Bej^£2(£',a;,f),i?(ti)]. After substitution of this commutator 
in (6.25) one gets 

-i5o,,,ME + io) [J2 ^^'-Si,e2 {E, w - oj', t)Ut + 5x,^^D^^Ji^\E)Ut) 
-ih,e,lx{E + a;) Dt^,Bo,e, {E, io - J, t)Ut + 5o,e,Dt^Ho\E)Ut) 

Ul' 

Prom this it follows that 

Bo,o(-B, t)Ut = UtBo,o{E, UJ, t) - i-fo{E + w) ^ D^,Bi,o{E, uj - uj' , t)Ut (6.26) 

Bifi{E, UJ, t)Ut = UtBi^o{E, UJ, t) - i-fi{E + uj) 

>^ {j2^-^'^o,o{E,u>- uj',t)Ut + Dt^^jL^\E)Ut) (6.27) 
After substitution of (6.27) in (6.26) one gets 
Bq^o{E, uj, t)Ut = -^Yl D^'DlL^"7o{E+(^hi {E+uj-uj')Bo,o{E, uj-uj'-uj", t)Ut 

ijj' Oj" 

+UtBofi{E, UJ, t) - i-fo{E + oj)J2 D^,UtBifi{E, uj - uj' , t) 

-^D^,D+_^jo{E + uj)^i{E + UJ- uj')iJio\E)Ut 

Now changing the summation index in the double sum, uj" = uj — uj' — uj" , using 

the definition of u>'{E\ and bringing the double sum to the left hand side of the 
equality, one obtains (6.21). The derivation of (6.22)-(6.24) can be done in a similar 
way. □ 

7. The Normally Ordered Equation 

In this section we will bring equation (6.20) to the normally ordered form. For this 
goal we will express terms like BE^^g^{E,u),t)Ut in the RHS of equation (6.20) in a 
form in which the annihilation operators are on the RHS of the evolution operator. 
This form is based on equations (6.21)-(6.24) of Lemma 1, which we will solve with 
respect to Bg^^g^{E,uj,t)Ut. 
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Suppose that for each e = 0, 1 the operators (1 + T^)ui.lj'{E) are invertible so 
that there exist the operators (1 + Te)'^]^„{E) with the propertie 

Uj" 

and that these operators are given by the convergent series 

oo 
n=l 

X E n,.SE)T^.,,..{E)...T^^_^^^.{E) (7.28) 



Detailed investigation of conditions under which this series converges will be done 
in a future paper. 
Then let us define 

rZ' (^) = E (1 + (^) 

R'jI'{E) = - E D^_^,{l+T,)zl^^{E)D+_^^^^{E + u;2) 

With these notations the normally ordered form of the equation (6.20) is given by 
the following theorem 

Theorem 7.2 The normally ordered form of equation (6.20) is 

dtUt=Y. J dE[Y,K\C'iE)J2l^^iE)B+,,{E,u;,t)UtB,,,,iE,u;',t) 

+ E [RIT iE)B+ {E, a;, t)Ut + Rl%'' {E)UtB,, {E, uj, t)) 

UJ 

+Rl]t{E) < ge,,PEe-^"'ge, > Ut\ (7.29) 

Proof. Using the inverse operators (1 + T^)Z\ji{E) in equations (6.21)-(6.24) one 
can express the products -Bei,£2 {E, t)Ut in terms of the products UtB^i^^^i^ {E, iv', t). 
For example: 

Bo,oiE, u, t)Ut = E(l + To)Z]u.. {E) \UtBQ,o{E, w', t) - i E D^'-^"lo{E + oj') 
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xUtBi^o{E,Lj",t)-J2Du."+^'D+, <go,PEL''go>ME + i^')7i{E-Lj")Ut 



and similarly for the other terms. Then after substitution of these expressions 
in (6.20) one obtains (7.29). □ 

It is known [1] that any normally ordered white noise equation is equivalent to a 
quantum stochastic differential equation. In particular equation (7.29) is equivalent 
to the quantum stochastic differential equation for the evolution operator: 



dUt=J2 J 



Y,K\:f{E)dN,{zi\cnE)) 



+Ro]t{E) < g,„PEe-^"'g,, > dt 



(7.30) 



Here we denote by (|/ >< (/|)a; the element of K^^ which belongs to K^^ sub- 
wen 

space. Moreover, 

Zl\':nE):^K^^^K^ 
and acts on an element X e Ki^n as 

zl\:^?{E)x = d^>,^„Y.fi,{E) < (Iff,, X PEge\h',x>^, (Iff,, X PEff.l)^ e 



8. Connection with Scattering Theory 

Let us show that the evolution operator in the LDL is directly related with the (1- 
particlc) T-opcrator describing the scattering of the system on one particle of the 
reservoir. Notice that the mean value of the evolution operator (2.1) with respect 
to the state (2.3) in the low density limit is equal to the vacuum mean value of the 
solution of the QSDE (7.30) 

lim (c/(VO) =< Ut >vac= e-" (8.31) 

where, since we average only over the reservoir degrees of freedom, the drift term T 
is an operator acting in the system Hilbert space Hs as 

r = - ^ / dERl%{E) < g,,PEe-^"^ge > 

£=0,1 

with RqIi{E) given at the beginning of section 7. 

Let us remind the definition of T-operator. For the interaction of scattering type 
the closed subspace of Hs ® r(Wi) generated by vectors of the form u (g) A+(/)$ 
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{u S 7^5, / G Til and $ G r(7ii) if' tlic vacuum vector) which is naturally isomorphic 
to Tis ® Til) is globally invariant under the time evolution operator exp [1(1/5 ®1 + 
1 (g) Hn + V)t]. Explicitly the restriction of the time evolution operator to this 
subspace is given by 

exp[i(ifs O 1 + 1 O ifi + Vi)t] e B{ns O Hi) 

with 

14 = ^ |5e >< (8.32) 

£=0,1 

The 1-particle M0ller wave operators arc defined by: 
f2±=s- lim exp[-i(iJs® l + l®-ffi + Vi)t]exp[i(iJs® l + l(g)iJi)t] (8.33) 

t— >itoo 

and the 1-particle T-operator is defined by: 

T = ViQ+ (8.34) 
From (8.32) and (8.33) it follows that 

f2± = s - lim U[^'' 

t— '■itoo 

where U^^'' is the solution of 

5t[/i^) = -iC/fVi(t), C/^'^ = i. 

Here 

Vr{t) = ^D,e'*- ® |5_tffo X S-M + ^i^+e-"- X S-M (8.35) 

We will show in Appendix that 

T = i I dEj2 R% {E)®\g,>< g,, \ Pe ■ (8.36) 

e,£'=0,l'' w 

Prom this formula it follows (see Appendix) that the drift term F is connected with 
the T-operator by the equality 

T = iJ2Pk {Tr-H,e-'^"'T) Pk (8.37) 

k 

This formula means that the drift term is the diagonal part (in the sense of the 
spectral projections of the system Hamiltonian) of the partial expectation of the 
T-operator in the 1-particle reservoir Gibbs state. In particular, under the RWA 
assumption (2.4) with wq = the expression for the drift term has the form 
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Formula (8.37) has important implications for the master equation which will be 
discussed in a future paper. 

9. The Langevin Equation 

Using the stochastic golden rule wc can find the Langevin equation, which is the 
limit of the Heisenberg evolution in interaction representation, of any observable 
X = Xs (S) lij of the system. The Langevin equation is the equation satisfied by 
the stochastic flow jt, defined by 

jt{X)=Xt:=U;XUt 

where Ut satisfies equation (7.30). To derive the Langevin equation we may apply 
the stochastic golden rule, described in previous sections. As a result we get in 
terms of the white noise operators 



Xt^Y. j '^^ 



+ T.{K,e,iE,co,t)U:%]^,'^{X)U, + U;&l^^^'^{X)U,B,,^,,{E,co,t)) 



(9.38) 



with the maps 

{X):= XRt,\% (E) + {E)X + 2 {E + u:)R+%\^ {E)XRl;% {E) 

S,UJ 

The equation (9.38) can be rewritten in terms of the stochastic difi'erentials as 



djt{X)=jtoY, J dE 



+ E(Q-T(^)'^^t^((l5^i X Pe9M + @l%'\E)dBt{{\g,, X PEge.lU) 

+jt o eo{X)dt (9.39) 



Here 



eo(X) := E / dEl'e'iE) [XRI%{E) + R+riE)X 
+2 E ^^le'{E + u:)Rly{E)XRil^{E) 



(9.40) 



is a quantum Markovian generator. The structure map Qo{X) has the standard 
form of the generator of a master equation [7] 



Qo{X) = ^{X) - -{*(!), X} + i[H, X] 
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where 



is a completely positive map and 
is selfadjoint. 



10. Conclusions 

An important problem in theory of open quantum systems is a rigorous derivation 
of a quantum Boltzmann equation. In the present paper the quantum model of the 
test particle interacting with the Bose gas has been considered. For this model we 
have presented a rigorous derivation of the quantum Langevin and the quantum 
master equation. 

We developed the stochastic limit method for the low density case. The pro- 
cedure of the deduction of a (unitary) evolution of the compound system in the 
limit of long time and small density of particles of the gas was developed. This 
procedure is being called the stochastic golden rule for the low density limit. The 
limiting evolution is directly expressed in terms of the 1-particle T-operator de- 
scribing scattering of the test particle on one particle of the reservoir. After that 
we obtain the quantum Langevin equation. This equation includes not only the 
system but also the reservoir dynamics. The master equation can be obtained af- 
ter averaging of this equation over the equilibrium state of the reservoir. We find 
the generator of this master equation which describes the reduced evolution of the 
test particle. We show that this generator has the Lindblad form of most general 
generators of completely positive semigroups. 

We considered the situation when the temperature of the reservoir is high enough 
so that no condensate is present. It is an interesting important problem to generalize 
the stochastic limit method to the case that Bose condensate is present. 



Appendix 

Let us first derive an explicit formula for the T-operator (8.36). 
The perturbation expansion for is 

^+ = y2{~^r / dh... dtMtn) . . . Viih) 

n=0 

This expansion induces the following expansion for the T-operator: 

oo 
n=l 
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/■OO rtn-l 

T„+i = (-i)" / dh... dtnViViitn) . . . Vi{h) 

Jo Jo 

and Vi{t) is given by (8.35). 

By direct calculations one can prove that 

T2n = iJ dE {T^^{E)\go >< Pe9o\ + T^l{E)\g^ >< PEgi\) 
T2n+i =iJdE (TO„Vi(i?)|ffo X PE9i\+Ti^+im9i >< Pe9o\) 

with 

T°0(£;) = (-1)" Yl • • • D.,r.-.D+,^.ME - ^2n-i) 

aJ,aJi...uJ2n-i 

X 7o(-B - U)2n-1 + W2n-2) • • • 7o(-E - U)2n-1 + ■ ■ -+^2)71 {E - W2n-1 + . . . + ^2 - ^l) 

T^^{E) = (-1)" • • • D^,n-.D.,n-^ME + 0J2n-l) 

u;,u;i ...aJ2Ti-i 

X 7i(i^ + W2n-l - W2„-2) • • •7l(^ + '^2n-l " • • ■ - <^2ho{E + (jJ2n-l - ■ ■ ■ - UI2 + <^l) 

T^n+l{E) = -*(-!)" E ^-^i • • • ^i„-i^-2„70(i? + UJ2n) 

X 7i(£ + a;2n - t^2n-l) • . .7o(-B + W2n - • • . + (^2)71 (-^ + W2n - . . . + ^2 - ^l) 

nn+l{E) = -i(-l)" ^^^-1 • • • D^.r.-M,ME - ^2n) 

X 7o(£; - L02n + W2n-l) • • •7i(-E^ - a;2n + • • • " t^2)7o(-B - a;2„ + . . . - ^2 + Wi) 
Let us show, for example, that 

00 

YR]jl{E) = Yn'n+.{E) (A.1) 

o) n=0 

In fact, 

YR'AiE) = -zY D^il + T,)-J^,o{E) = 

oj aJ,aJi 

00 

-iYD+6^„o-iYi-^r E D+T^u^^{E)T^^^^^{E)...Tl^,{E) 

ojfOJi n=l oj,oJi...u;n 
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One has 

U;,UJi...Ct'n 

(-1)" Y. D+T^u^^{E)T^^^^^{E)...T^^^_^^^^^^^i^ = 

W,a>l,W3...a)2n-l 

X 7o(-E + wi)7i(£-a;2) ■ . .7o(-E' + W2p-i)7i(-£^-'^2p) • . .jo{E + W2n-ihi{E - ui2n) 

Now let us make a change of summation index W2n-i W2n-i — W2n, then W2n-2 
0J2n-2 — W2n-i + W2n, •■• and finally u}\ ^ ijji — u)2 + ■ ■ ■ + t^2ri-i — W2n- After this 
change of summation indices we get for the RHS of the previous equality 

(-1)" D+D^,D+^...D^,^_,D+^^...D^,^_,D+^^^o{E+uj,-uj2 + ..--cv2n) 

X7i(£; - u;2 + - . . - C^2n) . . . 7o(-£^ + ^2n-l " u;2n)7l(^ " ^2n) = ^2"n+l(^) 

and from this (A.l) follows. One can prove in the similar way that 

oo oo 
uj n— n=l 

Therefore the T-operator is given by formula (8.36). 

Now let us prove the relation (8.37) beetwen the drift term and the T-operator. 
Since go and gi are mutually ortogonal one has 

iTrn,e-^"'T = - E / ^^E^-.o(i^) < 5e, Pi^e-'^^^s. > 

£=0,1"' OJ 

Then by expanding the operators (1 -l-Te)^^, in the series (7.28) one can show that 
for w 9^ 

E^fe^:',o^fe=0 
fe 

and 

E^fe^^o;o^fe = ^o;o 
fe 

Prom this (8.37) follows. 
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